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Use of RT-PCR and DNA Microgrrays to Characterize RNA 
Recovered by Non-lnvasSve Tape Harvesting of Normal and 
Inflamed Skin 

Rita. Wong * Vynga Tim* Vera Morhenn.t She^pin Hung-tf[ Bogi Andersen,§1f Elaine lto,fl[ Q. Wesley 
Hatfield, Ml anc| Nicholas R. Benson* * ™ ^ « r 

♦DermT-oti International Avanua of same, San Dlaao, California, MSA; iPWNan of DermBtoioQy, VA Med|a«| C#n^r r M Joll* VMam Dr, S«n Dtagpr California, 
uw!tKiSn^^ of California Irvine CallfprnH, USA; Mim** of m*QM Ch.mistr/, College of 

M^c|K IWwttM* * Own** Slolnforniatic., University of CeHIMnia, Irvine, California, USA 

Wo describe a noninvasive approach for recovering RNA from the surface of skin via a simple tape stripping 
procedure that permits a direct quantitative and qualitative assessment of petholofllc and physiologic Womarkers r 
Using seml-quantitatlvo RT^PCR wo show that tapeNharvflsM RNA i* comparable !n quality and utility to RNA 
recovered by biopsy, It is likely that taperharvestad RNA Is derived from epidermal cells residing close to (he 
surface and includes adne*a| structures and present data showing that tape and biopsy likely recover different cell 
populations. We report the successful amplification of tape-harvested RNA for hybridization to DNA microarrays. 
The*a experiments showed no significant gene expression level differences between replicate sites on a subject 
and minimal differences between a male and female subject We also compared the array generated RNA profiles 
between normal and 24 h 1% SJ-S-occluded skin and observed that SLS treatment resulted \n statistically 
significant changes in the expression levels of more than 1,700 genes, These date establish the utility of tape 
harvesting as a non-invasive method for capturing RNA from human skin and support the hypothesis that tape 
harvesting Is an efficient method for sampling the epidermis and Identifying select differentially regulated 
epidermal biomarkers. 

Key words: tape strip/dermatls/inflammatlon/DNA mlcroarray/quantitive RT-PCR/epldermie/expression profile 
J Invest Dermatol 123;159-167< 2004 



Contact dermatitis, a common skin reaction, involves 
several signaling pathways. Irritant contact dermatitis 
(ICD) predominantly Involves keratlnocyta activation (Freed- 
berg ef a/ T 2001), whereas Langerhans' cell presentation of 
antigen to Tcells in draining lymph nodes and recognition of 
the offending allergen in skin by memory Tcells control the 
Initiation and expression of allergic contact dermatitis (ACQ; 
Feghali and Wright, 1997) Clinically, both contact darmrtl- 
tides are characterized by pruritus, erythema* and edema. 
This commonality of the clinical signs and symptoms maKes 
distinguishing between ICP and ACQ difficult at the clinical 
level, particularly when symptoms are subtle, By contrast, 
at the molecular level, ICD and ACQ are believed to be 
characterized by unique rpRNA patterns, although the 
PMblished literature is conflicting (Hoefa.KKer ef a7, 1995; 
Filer era/, 1999; Morhenn eta/, 1999; Ryan and Gerberjck, 
1999; Utfgren ef a/, 2000; Cumberbatch et a/, 2PQ2), 
Documentation of simple and complex mRNA profiles la 
possible using RT-PCR and DNA microarray technologies. 
Using the technique of tape sfrlppinSr RNA can be 
harvested from both normal and Inflamed skin (Morhenn 
et a/, 1999) and by combining tape stripping and RNA 



Abbreviation: PPDE, posterior probability of differential expression; 
SLS, sodium lauryl sulfate 



profiling, it may be possible to non-lnvasively establish a 
diagnosis of ICD or ACD. 

In this study, we expand upon the work of Morhenn er si 
(1999) who used adhesive tape as a substitute for punch 
biopsy to non-invaslvely sample the epidermis, Morhenn 
eteil demonstrated that when a skin site was tape stripped 
20 or more times, it was possible to recover sufficient RNA 
from skin cells adherent to the tapes to detect and seml- 
quantlly specific mRNAs using the rlbpnudease protection 
assay. We demonstrate in this work that sufficient RNA can 
be recovered using sequential application of as few as four 
small tapes. 

In order to document the use of tape harvesting as an 
accurate and reliable sampling method w© performed a 
clinical trial in which occlusive patches containing erther 1 % 
SLS (irritant) or water (vehicle control) were applied to the 
mid-back of 10 subjects for 24 h, The sites were then 
clinically assessed and, along with normal control skin, 
surface cells were harvested with four applications of 
individual tapes and by shave biopsy. RNA was extracted 
from the tapes and biopsies and assayed semhquantita- 
tively for IHp, IL-8, GAPDH, and pectin mRNA using 
fluorescent, quantitative RT-PCR. The results showed 
consistent Increases in IU-1 p and IU-a mRNA In Inflamed 
skin relative to untreated skin. We further report the 
successful use of tape-harvested RNA to profile normal 
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and experimentally Inflamed sWn using DNA mlPrwws, 
This profile of SLS-|rrftatscl sKln Is the first step in the 
definition of RNA profiles designed to differentiate irritant 
from allergic skin reactions. 



Results 

Total RNA yields RNA was recovered from 27 of 30 sKip 
sites using four tapes as described In Materials and 
Methods, The amount of total RNA recovered was variable 
from site to site and subject to subject (data not shown). 
The average mass of RNA recovered from normal akin sites 
was 0.92 ng ( ± 0.35) with a range of 0 (two samples)-3 r a 
ng. The average mass of RNA recovered from water- 
occluded skin was 0.69 ng (±0,27) With a range of 0 (one 
samp|e)-2.7 ng- SLS inflamed skin produced the greatest 
average yield of RNA with an average of 1 85 ng ( ± 78) and 
a, range of 0,067-747 ng. 

Relative levels of housekeeping genes In tape-strip 
samples and biopsies We have chP&an aS markers of the 
Inflammatory process IL-1 p and IL-a mRNAs. We have 
accounted for differential recovery of total RNA mass In a, 
sample by normalizing these mRNAs to an Internal control, 
the (3-actln transcript. We then calibrated this normalized 
RNA ratio to the analogous ratio In a control sample (normal 
skin) using the comparative or A&G X method (described In 
Materials and Methods), In this study. IL-lp and mRNA 
are predicted to increase relative to pectin in response to 
SLS treatment, whereas the level of housekeeping mRNAs, 
such as p-actjn and GAPDH are assumed to remain 
constant. We have tested the assumption that house- 
keeping mRNAs such as GAPPH and p-actin are unchang- 
ing among different samples by measuring the relative ratio 
of these two mRNAs, 

If the relative ratio of two housekeeping genes is 
unchanged in different samples, the relative ratio of 



GAPDH/p-actin mRNA between two samples should be 
equal to 1 , The data In Table I reveal the fold-change of the 
QAPDH/p-actin mRNA ratio In SLS- and water-treated 
samples relative to normal skin and a|so water-treated skin. 
The tape sample data show that the average fold-change jn 
SLS samples was 0,55, whereas the average fold-change in 
water-treated samples was 1,14. Biopsy samples showed 
similar and minor changes with SLS-treated samples having 
an average 0.39 fold^cfiange, whereas water-treated sam- 
ples have an average o,53 fold-change. Individual subject 
data are shown in Online T&We S1 , Although there are some 
examples of statistically significant changes relative to the 
normal skin, these changes are not sufficient to explain the 
fold-changes jn IL-1 p and IL-a, which are reported In T&bla I, 
Thus, It Is likely that our housekeeping genes do change 
relative levels In response to SLS and water treatment, but 
the magnitude of those changes Is minor. Similar observa- 
tions and conclusions have been reported by Paludan and 
Thestrup-Pedersen (1Q92) and Grangsjo etGl (1996). 

IL-1p/p-actln mRNA ratios In SLS-lrrltatod and control 
skin Data in Table I reveal the average fold-change of the 
IL-ip/p-actin mRNA ratio In water-occluded and SLS- 
occluded skin relative to normal skin, In tape- and biopsy- 
harvested RNA samples. Individual subject data can be 
found In Online Table $2. Table I shows that on average the 
IL-1 p/p-actln mRNA ratio In biopsy samples of SLS-treated 
skin was elevated at least 11 -fold (>95% confidence 
Interval) compared with normal skin, Furthermore, the IL- 
1p/actin ratio was 4.42-fold elevated In SLS- relative to 
water-treated samples. In six of 1 0 subjects water occlusion 
produced significant (>95% confidence) Increases In the 
IL-1 p/actin ratio but this increase was typically 2-4-fold and 
was always smaller than the respective effect of SLS 
occlusion (Online Table S3). Thus SLS-occluslon produced 
the most consistent elevation of the M p/actln mRNA ratio 
but water-occlusion did effect a similar albeit smaller 
response. 



Table!. Summary of fold-clwnfle In GAPDH, IL-1p, and IL-8 mRNA relative to p-actln mRNA In SLS-treated and water-treated skin 







Average foid-chango by sampling method and treatment* 






Tape 


Biopsy 


mRNA 


Calibrator* 


Wat«r 


SLS 


Water 


SLS 


GAPDH 


Normal 


1,14 (O.e-2.14) 


0,55 (0,31-0.99) 


0,53 (0.43-0.63) 


Q,39 (Q,29-Q.53) 




Water 


1 


04A (0.36-Q.B5) 


1 


0.74 (0,55-1) 


IL-1p 


Normal 


NO 


>2.86 


>2.49 


>11 




Water 


1 


>1.S 


1 


4.42(1-42-13.76) 


IL-a 


Normal 


1,94 (0.93-5.90) 


10,32 (3.B9-27.34) 


>1,54 


> 52.45 




Water 


1 


6.31 (2-14.1) 


1 


34.08 (6,37-161.67) 



"Averaoa fold-change in mp to 1 o aubJeota (n m a for normal awn tape samples, n = e for wAter-trwated sKln tape sample* and n - 10 for a|| others). The 
Bverade fold-cha,nge Id C0|cu|at»cj from tho aubiBct-warago AC* values given In Online Tobies S3, S4. S7. Individual fo|d-chano* valuta are reported In 
Online Tables SI , S2, S6. Values preceded by > W lower llmft estimates n»cess|tated by th« fact that the applicable mRNA was not detectable in the 
control (calibration) samples. The > 95% confidence intarval for fpld-Cheng* ^ given In parflnthosls. 

"SLS and wsfsr samples are calibrated to normal skin aempfwi additionally, SLS samples are calibrated to the water sample. Calibration la described 
(n Materials and Methods. 
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Data in Table I also reveal the fold-change of l|_-ip/p^ 
actfn mRNA in tape-harvested samples of water^occlucjed 
and SLS-occluded skin relative to water-occl^dec) and 
normal skjn. IndMdqal subject data are shown In Online 
Tables S2 and S3. In most subjects, in p was undetectable 
in normal and water>treated, tape-harvested sKin sample, 
The average fold-change of the IL-1 p/actln ratio was 
estimated to be at least 2,66-fold In SLS-treated samples, 
Although the average fold-change for eight sublets was 
modest five of these subjects showed significant (>95% 
confidence) IL-1 (J/acttn ratio increases that were in quality 
tive agreement with the biopsy data (Online Table S3), 
Analysis of the remaining five samples was indeterminate 
because IL-1p mRNA was not detected in either the water 
or normal skin sample, 

IL-a/p-actin mRNA ratios In SLS-lrrlt9ted and control 
skin Table I reveals the fold-change of the |L-8/p-actin 
mRNA ratio In water-occluded and SLS-occluded s |<|n 
relative to normal skin. The data demonstrate that the |L- 
8/actln ratio was on average 52-fold elevated In SLS-treated 
skin samples relative to normal skin end 34-fold elevated 
relative to wafer-treated skin. Biopsy samples from nine of 
10 untreated skin sites had undetectable IL-B mRNA |evp|s 
and the sole normal skin biopsy sample with detectable II^b 
mRNA was close to the level of detection (Online Tables S4 
and £5). Thus IL-8 mRNA was generally not detectable In a 
biopsy of normal skin- 

Table I reveals that the IL-8/actin mRNA ratio was pn 
average 10-fold increased In tape-harvested samples from 
SLS-occ)uded sites (>QB% confidence interval) and 1-94- 
fold increased In water-treated samples {not significantly 
different than normal skin)- We conclude that the tape data 
are in good qualitative agreement with the biopsy data, with 
a majority of Inflamed sites revealing increases in 1^8 
mRNA. 

Purlng our analysis of IL-B mRNA in biopsy snd tape 
samples, we observed that although the fold-change was In 
qualitative agreement between the two methods, the 
primary ACj data were strikingly different (Online Table 
S4). As described In Materials and Methods, AC t is a 
measure of the ratio of two mRNAs in a sample. Online 
Table S4 shows that the AC tl | U . B for tape samples is very 
different than that of biopsy samples -for both SLS- and 
water-treated skin sltes-and this difference Is highly 
significant (both p-values < 0,005). The most striking 
example of this Is revealed by a comparison of the average 
AC liM a in tape-harvested samples compared with biopsy- 
harvested samples of water-occluded sites. For tape- 
harvested samples the average (ACt^-eWer is 1*64 r 
whereas for biopsy samples It is 9.22 (Online Table 54), 
The data suggest that In the call population harvested by 
tape, IL-8 mRNA Is much more abundant relative to p-actjn 
mRNA than in biopsy samples. TT>e average relative 
abundance can be estimated as equal to 2 _n,64 ^' 22 > or 
205, In this calculation, the biopsy sample acts as the 
"calibrator", thus the tape samples have, op average, a 
205-fold greater IL-6/actln mRNA ratio than biopsies. A 
similar calculation can be performed with the QLS sample 
data (Online Table 54), which also reveals a highly sig- 
nificant difference between tape and biopsy (AGuL-a^s 
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values (p< 0.005). We conclude from this data that tape 
harvesting recovers a distinctly different cell population than 
does biopsy, 

DNA mlcroarray analysis of RNA extracted from normal 
skin using tape The success In the above analysis of tape- 
harvested RNA from different skin sites suggested that this 
RNA might be amenable to amplification and hybridization 
to PNA microarrays. In order to assess the reproducibility 
and consistency of tape-harvested RNA samples for gene 
expression profiling experiments, we collected three sam- 
ples from the upper back of each of two healthy individuals, 
one male {sample C1, 02, and C3) and one female (sample 
A5, A6, and A9), Approximately 1 ng of total RNA was 
Isolated and the mRNA was amplified and blotin labeled 
using a MessageAmp aRNA kit (Amblon lnc„ Austin, Texas) 
as described In Materials and Methods, The resuming biotln- 
labeled aRNA from each sample was used for hybridization 
to an Affymetrlx HG-U133A GeneChlp, 

The results In Table |l ?how the differences observed 
when a matrix of pairwise gene expression comparisons 
between two GsneChlps was performed using Affymetrix 
Mlcroarray Suite software. Ttiese data show an average of 
only 12% variance among gene measurements, regardless 
of whether data from different sites on the same Individual 
or sites from different Individuals are compared, Further- 
more, comparing the data In quadrant three of Table II (A 
versus C) to the data in quadrants one (A versus A) and four 
(0 versus C) shows that about 1 b% of this variance Is due to 
either gender difference (A versus C) or inter-subject 
variation (A versus A or C versus C)> Thus, amazingly little 
variance is contributed by samples obtained from different 
sites or from different Individuals. 

To compare these data in a more quantitative manner, 
the three Affymetrlx GeneChlps each hybridized with targets 
from RNA samples obtained from individual A were 
compared with three QeneGhlps hybridized with targets 
from the three RNA samples obtained from individual C. 
These data were analyzed with a regularized t test (Long 
ef a/ ( 2001) implemented in the Cyber-T statistical program. 
This three-by-three comparison revealed 21, 790 probe sets 



Table II. Percentage of the measurement of gene expression 
reported as uncharged by the Affymetrlx MAS 5.0 software for 
each of nil possible palrwls* comparisons among GeneChlps 
(AS, AG, AO, C1, C2, nnd C3) hybridized with aRNA obtained 
from three different locations on the upper back of two 



subjects (A and C) 





GeneChlp/eubJeot |P (%) 


A5 


A6 


AO 


C1 


C2 


C3 


A5 


100 












A6 


ea.90 


100 










A9 


90.S0 


B6,1Q 


100 








C1 


60.80 


86.20 


67-40 


100 






C2 


B5.PD 


85,30 


83,10 


B9.60 


100 




C3 


88, GO 


86,00 


67.30 


86.00 


63JQ 


100 
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that exhibited gene expression levels above bapKground fpr 
all three sites from each subject. Of these genes 1117 (5%) 
were differentially expressed with p-value less than 0,0035, 
which based on the global false positive ana! negative levels 
of this dataset corresponds to a PPDE value of 0.95. Thus, 
56 of the 1117 Differentially expressed genes the* exceed 
this p-value threshold are expected to be false positives. 
The source of these inter-subject gene expression differ- 
ences remains to be determined; however, at least one of 
these differences is gender based. For example, the gene 
with the smallest p-va|ue and the highest PPDE value is the 
Y-linked rlbosomal protein S4 (PRS4Y). It Is likely that 
differences that are not gender based are a reflection of 
normal variation of gene expression between Individuals. 
These data are available at www.igb.ucl.edu, 

DNA mlcroarray analysis of normal versus water- 
occluded and SLS^occluded skin In a separate experi- 
ment, a total of nine RNA samples ranging from 1 to 10 ng 
were isolated by tape harvesting from three untreated, three 
water-occluded, and three SLS-occluded Sites of each Of 
three individuals. mBNA from each of the nine samples was 
amplified, blotin labeled and used for hybridisation to each 
of nine Affymetrix HG-U133A GeneChlps as shown In Fig 1 , 

Undated versus SLS-iraatetf samples A comparison of 
gene expression levels between threB untreated (A1, B1 P 
C1) samples and three SLSrtreated (A2, B2, C2) samples 
revealed 21,031 genes that exhibited expression levels 
above background for all samples, To assess the con^ 
fldence in global changes in gene expression, methods 
implemented In Cyber-T were used to determine the 
posterior PPDE of each gene based on experiment-wide 
global false positive and negative gene measurement levels 
as described by Hung ef eif (20Q2) and Baled and Hatfield 
(2002). When untreated versus S|_S-ocp|uded data are 
compared, the p-vatues for the differentially expressed 
genes are low and duster toward o. This |3 consistent with 
highly statistically significant differences among measure- 
ment levels of some genes. In fact 1771 genes are 
differentially expressed with a threshold of p- 0,003, WhlPh 
corresponds to an experiment-wide probability for differ- 
ential expression (PPDE value) equal to or greater than 0,90. 
These data are available at www.igb.uoj.edu. 

SLS- versus w$ter-treatod smples A comparison of gene 
expression levels between three SLS-treated (A2 F 02, C2) 
samples and three water-treated (A3, B3, C3) samples 
revealed 21,307 genes that exhibited expression levels 
above background for all samples. Based on a threshold of 
p =0,003, 1364 genes are differentially expressed with a 
PPDE value of 0,99, Of these, 1063 genes are also 
differentially expressed with a p-value of 0,003 and a PPPE 
value of 0.09 when SUS and untreated samples are 
compared. These data are available at www.igb-ucLedM. 

Untreated versus w$tert&Qt&cl samples A comparison pf 
gene expression levels between three untreated (A1, 01, 
Cl) samples and three water-treated (A3, B3, C3) samples 
revealed 21,164 genes that exhibited expression (eyels 
above background for all samples- This comparison 
revealed no statistically significant differential expression. 
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Nevertheless, based on our review of the genes assigned 
the lowest p-values. many of which are associated with 
Inflammation, we believe that the water treatment does lead 
to some changes In gene expression compared with 
untreated control skin. These data are available at www. 
Igb.ucLedu. 

For purposes of discussion, only the 100 genes dif- 
ferentially expressed With p-va|ues less than 1,4 x 10^ 
and PPDE values greater than 0.99 are discussed here, (n 
Fig 2 and are fully described in Online Table S6. An 
examination of these top 100 genes most significantly 
altered when the SLS-treated skin samples were compared 
with untreated skin samples revealed that, as expected, 
most of these genes carry out functions related to tissue 
Inflammation and injury (Fig 2; Online Table S6), These 
differentially expressed genes are proteinases, protease 
inhibitors, cytokines, chemokines, complement compo- 
nents, HLA factors, or receptors Involved In immune 
regulation. These associations with inflammation and injury 
responses for many of these mostly upregulated genes are 
documented in the literature (Online Table S6). These results 
demonstrate that the tape-stripping method described here 
harvests RNA suitable for complete genB expression pro- 
files of the skin that accurately reflect its pathological state. 



Discussion 

Recent advances in molecular medicine have made the 
possibility of molecular diagnosis a reality (Altman, 2001; 
Bertucci et al t 2Q01; Gallegue and Casellas, 20Q2; U*croix 
6t ah 2002; Whipple and Kuo, 2002; Satagopan and 



Individual A Individual ft In<tivtfuu) C 




IfcperitiMsnUi Alt M, CI wfttf A^B^CS 
Experiment 3 : A2, D2» C2 versus A3,ltt,C3 



Experimental 4*s|gn of SLS Irritation protocol, Details are described 
In the text. 
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complement And cqwhrmcut rccuprow 

poruxllMw niMrii iifracjaTtd hroWlrm 

lVwii«rtptlan ftcton 

OlJm*. 

* i4 ^ ( |i h w n If 1i 

F 'fl ure2 

Functional grouping of tho top 1QO differentially ™pr»wd genas 
between untmtic) and SLS-tmnt«d aKln with p^nlum lasa than 
1,4 x 1(T 1P and PPPE value* gn»«t»r than 0.99. Gemea tliat are 
unregulated In SI-3-traatad conditio^ gra mpreaantad with dark 
color and oerws th^t are dowryegulatad in sLs-traatad condltlpne 
represented by light gray. The number in parentheses |s the numby pf 
genae In each functional group that haa been previously documented to 
be related to the Inflammatory response. 



Panageas, 2003), Through the use of mlcroarrays and RNA 
profiling |t is becoming Increasingly clear that simple a,nd 
complex cell populations can be monitored or "profiW" 
with the intent of understanding the physiologic state pf 
those cells or tissues. This Information Is expected to lead 
to more accurate and possibly predictive diagnoses. We 
have shown here that the use of four small tape strips is an 
effective and basically non-invasive approach to capturing 
messenger RNA from the surface of skin and th*t this 
technique permits a direct quantitative and qualitative 
assessment of pathologic and physiologic blamarkers as 
a function of normal physiology. 

We have assayed seml-quantltatively the levels of IMP 
and IU-B mRNA relative to p-actin in normal, water and 5|_5<- 
occluded skin sites and shown that RNA from tape and 
biopsy samples produce qualitatively similar results. |n 
order to account for the possibility that changes In P-actin 
mRNA were responsible for observed changes In the 
interleukin/p-actln mRNA ratios (Suzuki of a/, 2000; Buatln. 
2002: TWcarfco et at, 2002) wo quantified tho levels of two 
housekeeping genes relative to each other. The resulting 
data (Table I) showed that although the QAPPH/p-actin 
mRNA ratio Is different In differently treated skin samples, 
the magnitude of this difference is not capable pf explaining 
the observed changes in Mp and iL-a mRNA levels 
Induced by S|_S treatment. This fact Is most clearly 
demonstrated by tape-harvest $nd biopsy data in Which 
IMP and IU-8 mRNA are virtually undetectable in control 
samples but easily detected In SUS-treated samples, an 
observation that cannot be explained by minor changes In 
p-actln mRNA levels. In addition, IL-1 p and IL-a mRNAs and 
proteins have been well characterized In Inflammation 
and are known to become elevated in response to S^S and 
other treatments (Paludan and Thestrup-Pedersen, 1992; 
Grangsjo ot 1996; Corslnl and Qalli, 1998; Tomlc^C^nic 



at al t 1998; Freedberg of a/, 2001; Perkins at af, 20Q1; 
Cumberbatch et fl/ r 2002; Coquette et a/ T 2003). Thus our 
data, from both tape and biopsy, are consistent with 
published observations, 

Biopsy and tape harvesting are not equivalent sampling 
methods, and therefore, should not be expected to yield 
Identical results. Tape harvest is restricted to the skin 
surface, and therefore, may preferentially recover vellus hair 
follicles and cells lining sebaceous, eccrine and sweat ducts 
as wall as cornoocytes (not predicted to contain RNA). Our 
method of using a single application of four Individual tapes 
dpes not result In glistening of normal skin end thus does 
not bare the viable epidermis. In contrast, a shave biopsy Is 
expected to include not only cells of the epidermis (primarily 
keratinocytss, melanocytes, and immune pells) but fibro- 
blasts from the upper dermis- The potential enrichment of 
surface epidermis and adnexal structures conveyed by 
our circular tape compared with a shave biopsy c^n be 
appreciated by considering that the surface arqa of a tape is 
284 mm 2 , whereas the surface area of a 2 x 2 mm shave 
biopsy is 4 mm 2 - 

The hypothesis that tape and biopsy recover different cell 
populations is supported by the AC ti|L .8 data previously 
discussed and presented in Online Table S4, Similarly, 
supportive data are presented In Online Tahl Q S7 for 
ACt.qAPDH values. Thus we propose that tape-harvested 
pells represent an enrichment of a subpopulatlon of cells 
found in a shave biopsy. Tills data imply that some 
differentially expressed blamarkers may be best detected 
in tape rather than bippsy-harvested epidermal samples. 

Identification of blomarKers diagnostic of clinical Irritation 
has been a long sought goal (Muller-pecker at at 1994; 
Boelsma ef a(, 1998; Muller-Decker etal, 1998; van Ruissen 
at a/, 1998; Komlne otst, 2001; Perkins et$l, 2001 J Boxman 
et a/, 2002; Perkins et a/. 2002; Coquette et a.1, 2003). We 
have used changes in IL-1 p and tlr8 mRNA w indicators of 
Irritation and shown that most but not all irritated sites 
display increased levels of these normalized mRNA markers 
(Table I and online Tables S2 and S5), Our data also show 
that tape-harvested and biopsy-recovered RNA are qualita- 
tively equal In their ability tP reveal an Irritant skin reaction, 
With respect to biopsy samples, it is evident that neither 
marker Is 100% efficient at diagnosing Irritation, a result 
observed for every biomarker proposed to be diagnostic of 
erythema and inflammation (Grangsjo et 4/, 1998; Muller- 
Decker et at 1998; Chung et s/ f 2001; Perkins et tf» 2QP1; 
Boxman at a/, 20Q2), The current limitation of the tape- 
harvest assay Is the inefficiency In detecting certain markers 
In samples with limiting amounts of RNA, a subject 
discussed below, But in comparison with the Sebutape 
assay (immunoassay of IL-8 protein; Perkins et a/ r 2001) for 
irritation, which has a sensitivity (Hoffrage et a/, 2000) of 
approximately 3096, mRNA biomarkers seem tc- possess 
superior potential- Our observation that water occlusion 
produced Increases of our biomarker ratios in some 
subjects has been reported by others {Grangsjo et a/, 
1996; Howie et a/, 1996; Perkins et */, 9001), 

It is reasonable to assume from the data presented here 
and in the literature that singular or limited multiplex assays 
will not be sufficiently predictive or diagnostic of irritation. 
Since our results show that the tape-stripping metfiod 
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harvests RNA suitable for DNA mlcroarray experiments, anc| 
that these gene expression profiles reflect the patMogJqaJ 
stats of human skin, however. It 3hogld be possible to identify 
a subset of ganes whose differential expression patterns can 
be correlated with different pathological states with a high 
degree of statistical accuracy. The fact that 1700 differentially 
expressed genes have beep identified with high statistical 
confidence sets the stage for tha creation of small custom 
PNA arrays designed to Identify patterns of gene expression 
diagnostic of irritant sKIn reactions, possibly diagnostic of 
different irritants and predictive of Irritant or toxic reactions- 
The next step along this path Is to identify the analogous gat 
of genes expressed during an allergic skin response, identify 
gepes unique to the Irritant or allergic response and combine 
them Into one DNA array, which could be used to determine If 
a mild reaction to a substance is Irritant or allergic In nature 
Such an array could also be used to test a variety of Irritants 
and allergens for unique profiles. 

Analysis of the top 100 genes differentially expressed in 
our SLS-freated samples shows fiat well over ha|f of thesa 
ganBS have been implicated jn injury and inflammation 
(Fig 2; Online Table S6), with most of these genes being 
upregulated. The DNA array data also confirm our RT-PCR 
data with IL-1 p being found In the top 100 most significantly 
differentially expressed (p-valqe less than 3.56x10 11 , 
probability of differential expression (PPDE) value >Q.99) 
genes, Although the 1^8 gene was not among the top 100 
most significantly differentially expressed genes, It m$ 
ranked number 395 In the normal versus SLS-treated 
comparisons (p^value = 6 x 10" e f PPDE = 0.a9) and ranked 
number 183 (p-value^5 x 10" 8 , PPDE-0.99) In the SL5* 
treated versus water-treated comparisons, The expression 
level of IH p was elevated 1 1 -fold and IL-B was elevated 0* 
fold in the 5LS-treated sample. Interestingly, many of the 
downregulated genes are hair keratins and keratin asso- 
ciated proteins selectively expressed in the hair during the 
anagen phase of the hair cycle. Either, the occlusive §LS 
treatment removes hair prior to th© tape stripping or the 
treatment blocks anagen in the hair follicles. 

We have shown that RNA cap be npn-invas|ve|y and 
productively recovered from the surface of the skin using toy r 
small tape strips. The number of tape strips can be reduced 
to two in conditions, where the surfacB of the skin has been 
disrupted, such as S|_S occlusion for 24 h (Nicholas Benson 
personal communication). Furthermore, the limitation of 
capturing small amounts of RNA from some akin sites can 
be effectively overcome by obtaining replicate control 
samples (data not shown) and by the appropriate choice of 
mflNA biomarker (discussed below). We have also presented 
data that the ACt^-a value, which is normally Msed to 
calculate a AAA value (calibrated fold-change), has revealed 
that tape and biopsy recover different cell populations. 

In this study, the quantity of PNA recovered from different 
individuals and skin sites was variable, with significantly 
more RNA being recovered from SLS-treated sites than 
normal skin sites. The large amount of RNA recovered from 
the SLS^Irritatsd sites, (s consistent with the known effects 
of sls, which effects invasion of Inflammatory cells end 
creates a weakened barrier facilitating the removal of the 
inflamed epidermis. We have since shown that recovery of 
RNA from normal skin is also a function of anatomical site 
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and similar sites vary between individuals With respect to 
RNA yield (data not shown). 

Although the variability of total RNA recovered doe3 not 
affect the results of relative gene quantitation, the recovery 
of very small amounts of RNA did affect our ability to fully 
analyze some samples, In this respect, the choice of 
blomarkere may be as important as the amount of RNA 
recovered from a site. For Instance, Onlina Table S7 shows 
that most tape-harvested samples could be assayed for 
pectin and GAPDH mRNAs and thus calibrated GAPDH/ 
actln ratios could be calculated (Online Table SI), Online 
Tables S2 and S5 reveal, however, that some of these same 
samples do not have calibrated M|Vactin or IL-a/actin 
mRNA ratios, with the IL-1 p assay being the most affected. 
The reason for this difference between biomarker assays 
lies in the relative abundance of tha specific mRNA. 
Because GAPDH mRNA is approximately equal in abun- 
dance to p-actln mRNA, all samples with detectable actln 
mRNA were successfully assayed for GAPDH- Likewise, we 
were highly successful at calibrating ILrB/actin mRNA ratios 
In water and SLS treated tape-harvested skin samples 
because |L-Q message Is relatively abundant In these 
samples. Thus the biomarKer mRNA that Is the most 
abundant relative to the normalizing mRNA will make the 
most efficient use of RNA mass. Therefore, candidate 
biomarker mRNAs should be chosen for best sensitivity, 
positive predictive value and high Native abundance when 
RT-PCR is to be used for detection and tape harvesting is to 
be the sampling method. 

We have chosen to demonstrate the utility of tape- 
harvested RNA for semi-quantitative RT-PCR and micro- 
array applications for several reasons. Both methods have 
particular advantages and are appropriate in different 
circumstances. The use of mlcroarrays Is an Invaluable tool 
for the discovery of diagnostic and prognostic biomarker 
candidates and may be essential for subcategortelng 
disease states, which may demand simultaneous assay of 
hundreds of blomarkers, The use of microarrays, however, 
[a expensive and technically laborious. Quantitative RT-PCR 
Is less expensive and less technically demanding and Is 
appropriate for studies where a limited number of known 
markers are being studied. 

In summary, our data show that the tape-stripping 
method non-invasjvely collects skin samples from normal 
and Inflamed skin that are suitable for RNA Isolation and 
gene expression profiling experiments. In the future, this 
method can be used to profile expression of a large number 
of genes In different skin conditions with the goal of 
designing custom arrays that allow molecular diagnoses 
of skin disorders. 



Materials and Methods 

cunlcal protocols The study protocols were reviewed and 
approved by an Independent IRP (BloMed IRP. San Piego, 
California)* all subjects efgnea: Informed consent and protocols 
adhered to the Declaration of Helsinki Guidelines, Ten. healthy 
women, agfls 21-55 were enrolled In the study, Regions of 
unblemished, normal appearing skin on the mid-hack were chosen 
for the application of two occlusive patches In the form of 
bandages approximately 4 cm x 6.5 cm, The bandages were 
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made using a clear, non^porous plastic hypoallargenic adhesive 
tape. In the center of this tape was a Weprll (non-woven cotton) 
patch measuring approximately 2 cm * 4.6 cm- One Wepril patch 
contained 0.P mL of1% aqueous sodium lauryl sulfate and the 
other contained 0.6 mL of sterile water as the veh(c|e control- 
Patches ware arranged such that the SIS patch was superior and 
directly adjacent to the water patch; the area of normal control skin 
was Inferior to and adjacent to the water patch* At 24 h pQSt 
application, the SLS and water patches were removed and the skin 
allowed to air dry tor 1 5 mln before scoring. The sites were scored 
py a trained technician uslnfl the sca|B provided Moy/. Patched 
sites were large enough that two areas could be tape harvested 
without overlap and room left for a shava Wopsy (^2x2 mm). 
SKfn sitae were tape stripped wlffl four tapes each] shave Mopsles 
were taken under local anesthetic flldocalne HC| 1% and eplnepfi^ 
Ine 1;100,0QQ; Abpott Laboratories, Abbott Port, Illinois), An area of 
normal skin was similarly tape harvested and shava biopsied, Tape3 
were stored In Individual appendorf tubes at -80 until extraction; 
biopsy samples were placed In buffer RUT and stored at -BO until 
extraction. Skin responses to each patch application were examinee- 
and graded under light supplied by a 100W Incandescent blue bulfr 
The following grading, scale was used: 0. no visible reaction; 1 , s||qrtf , 
pink, patchy erythema; 2, mild confluent, pink erythema; a. moderate 
erythema (definite redness) with edema: 4 T strong erythema (very 
Intense redness) with edema. In a second study, three subjects had 
three patches containing 1% SLS and three water patches applied 
to the m|d+ack far 24 h, Patches were removed, scored and tape 
stripped as above. In a third study, two Individuals were tape 
stripped on normal skin on the upper back at three adjacent sites as 
above, RfslA harvested In these lest two studies was used In the DNA 
microarray experiments described pelow. 

Materials and reagent* Adhesive tape was purchased from 
Adheslves Research (Glen Bock, Pennsylvania) In bulk rolls. These 
rolls ware custom fabricated into small circular discs, 17 mm In 
diameter, py Diagnostic Laminations Engineering {Oceanslfle, 
California), TotaJ apleon RIMA was pgrehased from Amblon, 
"RNeasy" RIMA extraction kit and Senslscrlpt Reverse Transcrip- 
tase kit were purchased from Qlagen (Valencia, California). FOR 
primers and probes OaqMan Pra-Developed Assay Reagents) §nd 
TaqMan Universal Master Mix, which Included all buffers and 
enzymes necessary far the amplification and fluorescent detection 
of specific cDNAs, were purchased from Applied BJoaysterns 
(Foster City, California). Total mRNA was amplified using the 
MessageAmp aRNA kit purchased from Amhlon Inc. (Au#|n, 
Texas), Human Genome U133A DNA chips were purchased from 
Affymetrlx Inc. (Santa Clara, California), 

Isolation of RNA The ANA within skin cells adherent to the four 
tapes gsed to harvest a site was pqo|ed by simultaneously 
extracting the tapes. In a volume of puffer RUT (supplied with 
RNeasy klt)r Extraction was performed using the manufacturer's 
directions and Included a Proteinase K digestion, sonlcatlon Of 
tapes and M on-column" DNase I digestion. RNA was e^ed In 100 
um of sterllo, RNase-free water. Extraction of biopsies was 
performed with the same kit according to the manufacturer's 
Instructions, 

Quantitative RT-PCRi Ten ^Lof RNA was reverse transcribed (RT) 
Into cDNA with the Senslscrlpt Reverse Transcriptase k|t ajyd 
random hexamers In a final volume of 20 |*L according to the 
manufacturer'3 directions. The reaction was diluted 5-fold With 
sterile, niiclease-free water (Amb)on) for use In the subsequent 
amplification/detection reaction, For each specific mRNA decep- 
tion, three replicate RT + reactions and one fTT (no PTi negative 
control) reaction were performed. Two ampllflcatlon/deteqtlpn 
reactions were done on each RT f reaction to yield a total of e|x 
Independent determinations of *e threshold value (C t ; diseased 
below). All RT reactions were amplified using two replicates and 
were negative (data not shown). 
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Quantitation of RNA mass recovered w(tfi adhealve tape and 
biopsy The amount of RNA recovered by tape Is too small On most 
samples) to detect by UV. We have also found that contaminants In 
the adhesive co-purify with the RNA and Interfere with UV and 
fluorornetric detection. We therefore estimated the RNA mass 
recovered from tapes py using quantitative RT-pCR with reference 
to a standard curve (C^in versus |og[RNAj; Applied Blosyatems, 
2001) created from commercially purchased human spleen total 
RNA. Spleen RNA was treated with DNase I and purified with the 
Qlagen RNeasy kit following the manufacturer's Instructions. 
Purified standard RNA was quantified spectroscopic^ |y using 
0,Di 260. The standard curve was constructed using four 
concentrations of RNA from 0,01 to ^ ng par mL Each RNA 
standard was reverse transcribed In triplicate ana) each RT reaction 
amplified once to yield three replicates per standard concentration. 
Amplification and detection of unknowns was accomplished as 
described below using p-actln mRNA as the quantified marker. 
Experimental samples were reverse transcribed In triplicate and 
each RT reaction amplified In duplicate to yield a total of six 
replicates, The average Ct, -cfln of these s|x replicates was used to 
calculate the concentration of PNA In the unknown with reference 
to the standard curve. The accuracy of this method relies on the 
relative amount of p-actln mRNA to total RNA In the epidermis 
being similar to that In human spleen. |f the relative amognt of p- 
actln mRNA to total RNA la different between the two tissues then 
our mass data will be similarly affected. Therefore, we describe all 
tape-harvested RNA mass calculations as estimates to reflect this 
uncertainty. RNA recovered from biopsies was quantified fluor- 
omeirlcally with the RlboGreen RNA Quantitation Reagent (Mole- 
cular Propes, Eugene, Oregon). 

Ampliflcatlpn and detection of specific mRNA Specific mRNAs 
were converted to cPNA as described apove. Specific pDNAa 
were seml-quantlfled using gope-speclflc primer/probes (5'-nuc|e- 
ase assay) and fluorescence detection- Amplification and detection 
assays were performed using TaqMan Pre-Deve|oped Assay 
Reagents (PDARi Applied Blosystems) on an Applied Blosysterns 
7000 Sequence Detection Syeterrip (L-1p mRNA was serriL 
quantised relative to p-actln In a multiplex assay where p-actin 
and IL-1P mRNA were detected simultaneously In the same Jupe, 
Similarly, IL-6 was seml-qMant&ated by multiplex PGR with p-actln- 
The data from these experiments are summarized In Table I and 
Individual subject data are reported In Online Tables S1-S5 and 57, 
All data were confirmed by slngleple* PGR assays where lL-fi, IL- 
1p and p-actln were separately assayed in Individual tubes (data 
not shown). GAPPH was seml-quantitated relative to p-actln using 
a slngleplex format. Thermal cycling conditions were; prior to 
cycling, 2 mln at B0°C f than 1 0 mln at 95*0; then 40 cycles at 95 fl C 
for 1 5 s and 6Q P C for 60 s. Threshold detection was set at 0,2 for all 
assays. 

Seml-quantltatlon of mRNA using the AAC, method and 
calculation of fold-lncreeae relative to a control sample In this 
work, we use the comparative or AAC ( method of calculating 
relative gene expression levels between two samples. In the AAC ( 
method the levels of IL-1P, M, and GAPDH mRNAs are assayed 
quantitatively and normalized to p-actln mRNA to create a ratio of 
the RNA of interest (I.e., M P, IL-0, QAPDH) to p-actln mRNA for 
each RNA sample. This ratio Is than further normalized to a control 
sample (''calibration"; Applied BiosyefemB, 2D01), This relative 
level of the RNA of Interest to p-actln |n a sample Is given hy the 
equation \i2"^ where *V Is an unknown constant related to the 
PGR conditions a,nd the fluorescent probe; because k Is unknown, 
the absolute ratio In a sample Is also unknown. ACt Is equal to the 
difference between the experimentally derived C t values for the 
mRNA of Interest and p-act|n mRNA (for Instance, A^oapdh" 
C t GAPPHr c i.»«n for GAPDH normalized to aetln). Although an 
absolute mRNA ratio cannot be determined for a sample with only 
C, data, the relative ratio petween two samples can be determined, 
This Is Illustrated below using, as an example, QAPDH as the 
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mRNA of Interest, pectin mRNA ae the natm^lzlno marker and 
calibration of the SLS sample to the normal sWn sample, TN ratio 
of gapph to actio mBNA In an SLS-treated sample Is given byi 

( mPNAoAPDn N _ j^c^ - ^"(Ot^wh-^Jm 

Likewise the ratio of GAPPH mflNA to actln mRNA In the n°rrnal 
(control) skin sample Is given Py 

AtiRNAgapdh \ 
VmRNAp^ftfitm,/^, 

By division of the first equation |ay the second, simplification, apd 
realization that under Identical PGR conditions the ponstants K ftre 
equal, It can be shown that I 

/ mRNAoApDHN Z/ mRNAqAppH N 
^ mPNAapiin / eu3 / V mRNAactin )^mv^\ 

where AACt,gu3 = AC^u* - AC, inW q, and ACys are defined aboU 
The C t values are the experimentally determined number of PpR 
cyc|Bs required to achieve a threshold fluorescence fatatlstlca||y 
significant Increase In fluorescence over background) for the Rf*|A 
of interest and p-actln mRNA (Gibson ef «/, 1996; He|d &t a/, 1998), 
Key to the comparative method and accuracy of Equation (3) |s 
the amplification efficiency qf the mRNA Mng quantified, the 
efficiency of pep must be very similar if not Identical for the 
ampllcona being assayed, In order to reassure ourselves that thlp la 
true In this study we Initially tested all the prlmer/prat^J sets at four 
different spleen RNA concentrations » In a plot of C t versus log[RNA 
mass] the slope of the line Is a. function of the efficiency] of 
amplification (Applied Plosystems, 2001). In these experiments, 
the slopes were essentially Identical and Equation (3) accurately 
portrays the relative ratios of our mRNAs (data not shown), 

RT-PCR data analysis Experimental da,ta are reported as the 
number of cycles (Cd required to achieve a threshold flucrascerp. 
Each reported C t Is the mean of a|x replicate, mea^remsnte. 
Calibrated fold-change calculations are made using Equation ft), A 
calibrated mRNA ratio is considered to be significantly different 
(>b&% confidence IntervaO than the control If the two standard 
deviation range of fold-change does not Include the value of 1, 
which la the defined va|ue of the calibrator (AAC t for the calibrator 
Is equal to 0» with a consequent fold-change relative to Itself eqw\ 
to 2\ The significance of AC» values (On line Tables S3, S4 t and^A 
was determined by applying a two-sided, pajred t teat. In some 
normal skin samples, IL-9 or IU1 p mRNA Is undetectable. In the.ee 
cases, we report an estimated lower limit of fold-change basecj on 
a lower limit of defection (C,=37). It Is our experience thfltittie 
practlcal limit of detection for reaHlme PGR le 37 cycles, When the 
threshold number of pep cycles (0^ extends beyond 37, C t v^MPs 
become highly variable or fluorescence does net achieve a 
threshold velqe, With sgch data, we assign a value of Ctf 37 
and consider the mRNA undetectable. To Illustrate a calculation of 
the lower limit of fold-change at >GS% confidence we use <Jata 
from subject 7, Subject 7'e tape sample of an SLS-treated site has 
a (ACt, |L-fl)eua - - 0.68 4 O,0B, The control AC t ^ (C^ iu^-A,«Wrwnw! 
= (Cnu^-2B^): because the |L-8 mRNA was undetectable we 
assign it a value of 37. hence AC )L * |n the control sample |s 37- 
29,6^7,2 (±0,13), the value reported In Online Table 64.,We 
can now calculate the calibrated fold-change of IL-8 in the SL5 
sample as 2-< MCt >=2 H -°' flp - 7 ' s » = 236, The lower limit of I the 
>95S4 confidence Interval, however, |s given as pr'^" 1 ' 3 ^ 
2 -<-7.*a ^2 * o.i« = ^ ^ va|ue sported in online Table S4 (P,1 8 
Is the error In AACt ((0.CN3) 2 + (0,1 3) 2 ) 1 ^). I 

I 

T7 linear RNA amplification mRNA wa3 amplified and hlptln 
labeled using a MessagsAmp aPNA kit purchased from Amnion 
Inc, according to the manufacturer's Instructions, Typical yields of 
aRNA obtained from two rounds of amplification ranged from 3Q to 
60 fig, 
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Hybridization "f blotlnylatad mRNA target* to Affymetrlx 
GeneChlp*. atalnlrm, data acquisition and data analysis Hy» 

brldteatlon and staining wer« performed according tp the manu- 
facturer's Instructions; a detailed protocol Is available at the 
manufacturer's website, 

DNA mlcroarray data ana|yil* The PIMA mlcroarray data re- 
ported here were globally normalized and pepkground subtracted 
with the Asymetrix MAS 5,0 software package and modeled with 
the dChlp program of LI and Wong K<M11, Statistical analyses were 
performed with fie Cyber-Tatatlstloal program aWllaMe on line at 
www,|gp,uc|,edu, The posterior probabilities of differential expres- 
sion (PPDE) values (ranfllpg from 0 to 1) for each gene were based 
on estimates of experlmen>wide false positive and negative levels 
calculated by Cyper-T according te the methods of Allison and 
Qadbury (2002) as described py (Hung ef a/ (200£() and Ba'dl and 
Hatfield (2GQ2), 
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